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Abstract. We present first results of the solutions of the Yukawa model as a Quantum Field Theory (QFT) 
solved non perturbatively with the help of lattice calculations. In particular we will focus on the possibility 
of binding two nucleons in the QFT, compared to the non relativistic result. 

PACS. 13.75.Cs Nucleon-nucleon interactions - ll.10.-z Field theory 



1 Introduction 

Nucleon-nucleon (TV TV) interaction is probably one of the 
most studied problems in theoretical physics. From meson 
exchange models |2E] till effective chiral lagrangians 3 , 
the effort of physicists has been towards the development 
of suitable NN potentials that, once included in a Lipmann- 
Schwinger (LS) equation, would provide the nuclear bind- 
ing energies and scattering properties. Using Green's Func- 
tion Montecarlo one can even compute the nuclear spec- 
trum of nuclei up to ~12 nucleons [2]. 

Most potential models are inspired by an underlying 
Quantum Field Theory (QFT), from which only a very 
particular kind of diagrams are taken into account when 
solving the LS equations: in practice the resolution is cur- 
rently available only in the ladder approximation. 

All crossed-ladder graphs were summed up for the Wick- 
Cutkosky (WC) model [3], and the resulting binding ener- 
gies are much bigger than those obtained within the ladder 
approximation. This strong bias is one of the most impor- 
tant motivations for the present work. As the WC model 
is not consistent as a field theory [5], we will study the 
simplest renormalizable QFT involving fermions, where 
one species of fermions interact with a scalar meson via a 
Yukawa coupling. 

The interest of this approach is manifold. On one hand 
it allows a comparison with the results of the ladder ap- 
proximation in different relativistic and non relativistic 
equations. On the other hand, and including other cou- 
plings, it could provide a relativistic description of nuclear 
ground states in terms of the traditional degrees of free- 
dom - mesons and nucleons - with no other restriction 
than those arising from the structureless character they 
are assumed to have. 
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2 The model 

We consider a system of two identical fermions (ip) in- 
teracting through the exchange of a scalar meson (</>) de- 
scribed by the lagrangian density, 

C = 4>Di> + C K g{4>) + gti$M> > (!) 

where D = 7^9^ — Mq is the Dirac operator with a bare 
fcrmion mass Mo, Ckg is a Klein-Gordon lagrangian for 
the scalar field. In the NR limit (|TJ gives rise to the po- 
tential 

„2 -m a r 

V{r) = - 9 f- , (2) 

47r r 

where m s is the meson mass. The NR model does depends 
on a unique parameter, G = f^r^-, and the first bound 
state appears for G w 1.68. The existence of this unique 
scaling parameter G can be easily shown in Schrodinger 
equation but it is no longer true for the relativistic case 
or the QFT. 

In order to study the bound states, one needs to take 
into account contributions to all orders in the coupling. 
Therefore, a perturbative approach is not suitable. In- 
stead, a genuinely non perturbative tool will be used, lat- 
tice field theory, developed in the context of QCD. The 
lattice Yukawa model is solved in a Euclidean space-time 
where vacuum expectation values are computed in the 
Feynman path integral approach. For the dressed nucleon 
propagator one has, for instance, 

G a ^x,y) = {0\r(x)^(y)\0)= (3) 

where the euclidean action acts as a probability distribu- 
tion, allowing for a Montecarlo integration. 
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We have chosen the following discretization of scalar 0,4 
fields: 



+ a 2 m 2 s )4>l-2j2<l>x+ 



(4) 



0,2 



and for fermion ones: 



xy 



where D xy is the Wilson-Dirac operator: 
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D xy = (1 + 9l4>x) S x , y - 

- [C 1 - T/i^x+Zi,, + (1+7m)^-A,»1 • ( 5 ) _0 ' 4 
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Fig. 1. Spectrum of Wilson-Dirac operator for a small lattice 
in which the hopping parameter, K = 1/(8 + 2aA/ ) and and several values of the coupling. 
ql = 2k<7o have been introduced. 

Fermion fields - being Grassmann variables - have to 
be integrated out in an algebraic way, resulting into: 



G at3 (x,y) = 1^ MD- l f y det(D)e- s ^ 



(6) 



This calculation is rather demanding in computing time 
due to the determinant. The task is considerably simpli- 
fied in the "quenched" approximation, which consists in 
neglecting all virtual nucleon-antinucleon pairs originated 
from the meson field (f> —> ipip. Because of the heaviness of 
the nucleoli, this appears as a good approximation for the 
problem at hand and has been adopted all along this work. 
Note that this is not a priori justified for QCD, where 
quarks are very light. Nevertheless, the quenched approx- 
imation gives there qualitatively good results. Mathemat- 
ically it is equivalent to set det(D) = 1. The main nu- 
merical task in calculating ||JjJl is the inversion of Dirac 
operator D xy 

In the quenched approximation, and in absence of me- 
son self-interaction terms, the meson fields are free, and <f> 
field configurations can be independently generated by a 
gaussian probability distribution in momentum space. 

3 Spectrum of Dirac operator 

The spectrum of the Dirac operator (JjJJ in the free case 
lies in a circle centered in A = (1, 0) and with radius 8k. In 
QCD when the interaction is tuned up, the eigenvalues are 
modified, but its real part is always bounded from below. 
In the Yukawa model, on the contrary, the coupling term 
plays the role of a mass: as the coupling constant grows, 
the spectrum spreads out in real part and some eigenvalues 
go to the negative real part half-plane (figure . 

Negative real part eigenvalues spoil the convergence of 
most iterative algorithms, but is not a fundamental prob- 
lem. For larger values of the coupling and large lattices, 
nevertheless, the probability of having one -or more- eigen- 
values very small grows dramatically. A simplified but sig- 
nificant picture can help to estimate the appearance of 




Fig. 2. Average number of negative eigenvalues for several 
lattices with Lams — 5 as a function of gh- 

those small eigenvalues. The diagonal terms in 10) have the 
form 1 + gL^x, that will be zero as soon as <p x = —1/gi,. 
The values of 4> x are distributed according to (TJJ, as a 
gaussian of width 
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2£(l-cos(*„)). (7) 



The probability of having one negative eigenvalue is plot- 
ted in figure [3 for different lattice sizes, showing how neg- 
ative eigenvalues appear for ~ 0.5. The non diagonal 
terms in |J5J modify this picture, and in practice there are 
small eigenvalues for g^ > 0.8, hindering the numerical 
solution of the linear system. This implies that there is a 
maximum value of the coupling constant that can be used 
in this model. The problem could perhaps be solved in 
the unquenched CcLSG, clS the fermionic determinant would 
eliminate the configurations with very small det(-D). 
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Fig. 3. One and two body effective masses as a function of 
the euclidean time, for a 20 3 x 80 lattice with go ~ 1.6 and 
am a = 0.200. In blue two times fermion mass, to compare 
with mass of two fermion + state. 

4 One and two-body masses 

One-body mass is computed from the time dependence of 
euclidean correlators: 

C 1 {t)=Y J Tr[G{x,t)]~e- M - t , (8) 

X 

for large values of t, determining the fermion renormal- 
ized mass, Mi. Preliminary results on one body masses 
for both scalar and pseudoscalar coupling were already 
presented in [Hj. 

Two body masses, M%, are obtained in a similar way, 
from the time evolution of the propagator of an operator 
J(x) — r a pip a (x)ipi3(x) creating a nucleon pair, 

C 2 (t) =^Tr(J(x,t)jt(0,Q)) ~ e- M2t , (9) 

X 

that for large values of t projects on the lowest energy 
state with the quantum numbers of the operator J(x). The 
matrix F determines the spin and parity of the state, being 
r = 1727075 for J* = + (ground state) and _T = 7270 for 

,r = rr. 

The exponential behavior is reached only at large val- 
ues of t. It is useful to define an effective mass as: 

that tends for large t to the mass of the state and helps 
to find the adequate fitting window. Some results can be 
found in figure for one and two body masses. 

In the figure twice fermion mass is plotted for com- 
parison. With these parameters twice fermion mass is not 
distinguishable of two fermion mass. This is the common 
picture for the whole set of parameters tested, no signal of 
the existence of a bound state is found below the critical 
value of the coupling co nstant . 

1 There might exist bound states for very light mesons, near 
the Coulomb limit, but this regime is difficult to reach on the 
lattice due to the hierarchy of scales appearing. 



5 Discussion 

Renormalization effects have been analised for one body 
masses, where perturbation theory works, and the renor- 
malization issues concerning the coupling constant have 
also been discussed [£]■ 

The existence of a maximum value of the coupling in a 
QFT treatment of the Yukawa model has been established. 
This critical value is smaller than the one needed to form a 
bound state in the NR limit, and no signal of such a bound 
state for lower couplings has been observed. This limit on 
the coupling is characteristic of the QFT, different from 
the potential approach where the coupling usually take 
values G>1. 

The meaning of this result needs to be clarified. It may 
be related to the quenched approximation, or the fact that 
we neglect meson self interaction. But then it should be 
noted that the same approximations are performed in the 
non relativistic (Schrodinger) treatment. For a given value 
of the lattice spacing, there exists other ways to discretize 
the nucleon-meson interaction which don't have these zero 
modes. This has to be further studied and it is not clear if 
it allows to reach larger renormalized coupling constants 
and particularly to reach the bound regime. We are not 
yet in a position to decide if the bound on the coupling 
constant we encounter is a lattice artefact or if it really 
casts a doubt on the Yukawa theory itself. 

It is known that the Yukawa theory is infrared free 
and, as such, encounters the "triviality problem" i.e. that 
the ultraviolet cut-off can not be driven to infinity without 
the theory becoming trivial. It means that this can only be 
an effective theory with a physical ultraviolet cut-off. The 
problem encountered from the difficulty to invert the Dirac 
operator seems also to put a limitation on the continuum 
limit. It is not clear whether both problems are related or 
just happen both to hinder the continuum limit. 

Whether we manage or not to overcome the difficulty 
of reaching the domain where bound states appear, the 
connection between the QFT treatment and the Schrodinger 
approach can still be performed by an estimate of the scat- 
tering parameters which can be computed thanks to the 
method proposed by Luscher [7|, and have recently been 
reexamined in 001]. 
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